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Abstract 

We propose a new scalar-tensor model which induces significant deviation from general 
relativity inside dense objects like neutron stars, while passing solar-system and terrestrial 
experiments, extending a model proposed by Damour and Esposito-Farese. Unlike their 
model, we employ a massive scalar field dubbed asymmetron so that it not only realizes proper 
cosmic evolution but also can account for the cold dark matter. In our model, asymmetron 
undergoes spontaneous scalarization inside dense objects, which results in reduction of the 
gravitational constant by a factor of order unity. This suggests that observational tests of 
constancy of the gravitational constant in high density phase are the effective ways to look 
into the asymmetron model. 


^This work was initiated and largely advanced while the author was staying at RESCEU as a visiting professor. 



1 Introduction 


General relativity(GR), which describes gravity in terms of a massless spin 2 field, has been 
tested and passed all the precision experimental tests such as the solar system and the terrestrial 
experiments [T]. No observations which clearly contradict with predictions of GR have been 
found. This does not guarantee that GR remains valid under extreme physical conditions beyond 
the present experimental limits. Indeed, it is known that GR cannot tell what happens at the 
center of black holes and at the very beginning of the Universe if energy conditions are to 
be satisfied [2]. In addition to this, the explanation of the accelerating expansion of the late 
time universe may require the modification of GR on very large scales. Motivated by these 
considerations, GR may be viewed as an effective theory which is valid only in some domain of 
space spanned by physical parameters such as length, energy and density, although the boundary 
of such domain is not yet well defined. Many possibilities have been proposed in the literature 
in various contexts (see references in m)- 

In the near future, direct detection of gravitational waves will become possible by using 
the laser interferometers such as advanced Laser Interferometric Gravitationalwave Observatory 
(aLIGO) [7], advanced Virgo (aVirgo) [8] and KAGRA [9]. Target gravitational waves originate 
from the vicinity of the compact objects such as neutron stars where the matter density is much 
larger than any other places in the Universe. Observation and analysis of such gravitational waves 
should enable us to probe the laws of gravity in such domain that has heretofore unexplored. 

Scalar-tensor(ST) theories are well-studied and natural alternatives to GR [ToHl^ . Obser¬ 
vations of gravitational waves enable us to probe ST theories in high density and strong gravity 
regime. Interesting targets relevant to gravitational wave observations are a class of ST theories 
which mimics GR in low density (or weak gravity) regime but significantly deviates from GR in 
high density(or strong gravity) regime |14lll5j . One natural way to construct such a model is to 
introduce interaction between the standard model particles and the scalar field by the conformal 
factor so that the effective potential for the scalar field depends on the matter density. If the 
system is static, the scalar field takes a value that minimizes the energy of the system. This 
expectation value depends on the matter density and controls the interaction strength between 
the standard model particles m- Then, it is possible that the expectation value vanishes if 
matter density is low and the spontaneous scalarization occurs if matter density exceeds a crit¬ 
ical value. In such a case, modification of GR occurs only in the high density region exceeding 
the critical density. This mechanism is completely opposite to the symmetron model proposed 
in m in which scalarization occurs only when the matter density becomes smaller than the 
critical density. In that model, matter density inside the solar system is supposed to be larger 
than the critical density and GR is recovered but deviation from GR appears on cosmological 
environment due to the low background density. For this reason, we call the scalar field that 
acquires non-vanishing expectation value only in high density environments asymmetron. 

In this context, there is an interesting scalar-tensor theory proposed by Damour and Esposito- 
Farese (DEF) |18lll9j in which significant deviation from GR occurs only in the vicinity and the 
inside of neutron stars and safely passes the solar system experiments. In the DEF model, the 
scalar held in high matter density region becomes tachyonic due to a particular form of the 
conformal coupling with the standard model particles (see left hgure of Fig. [T]). As a result, 
the scalar held takes a large non-vanishing value inside the neutron star and approaches a non¬ 
vanishing but much smaller value at distance far away from the star. The value at inhnity 
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Figure 1: Density dependent effective potential of the scalar field for the Damour-Esposito-Farese 
model (left figure) and the asymmetron model we consider in this paper (right figure). Dotted 
curve in each figure represents the value of 4> that (p would take inside the compact object with 
density p. 


is hxed to match the cosmological value just as in the case of the Fierz-Jordan-Brans-Dicke 
theory |2r)1123j and this value must be small enough to satisfy the solar system and terrestrial 
observational constraints. Since the magnitude of the scalar held controls the amount of deviation 
from GR, signihcant deviation from GR occurs only in the inside or the vicinity of the neutron 
stars. Because of this, the structure of the neutron stars differs from that under GR and this 
suggests that studying neutron stars and deriving observational consequences is the most effective 
way to test DEF model, as has been pursued in the literature [21H33]. 

However, it is known that the DEF model faces difficulty of embedding it in the cosmic 
history [28l|3lll35] . During inhationary and matter-dominated epochs, the coupling between the 
scalar held and the matter held forces the scalar held to take non-vanishing value and the law 
of gravity in the present Universe deviates from GR to an extent incompatible with the existing 
constraints. Our main motivation of this paper is to extend the DEF model to incorporate it in 
the cosmological context. We achieve this by dropping the two restrictions imposed in the DEF 
model. The hrst is the mass of the scalar held and the second is the energy scale appearing in the 
conformal factor. In the DEF model, the scalar held is assumed to be exactly massless and the 
energy scale in the conformal factor is taken to be around the Planck mass. We do not impose 
these conditions and assume that the scalar held is massive and the energy scale in the conformal 
factor differs from the Planck scale. Because of these assumptions, the effective potential of the 
asymmetron has a global minimum for any value of the matter density p, whereas the effective 
potential of the original DEF model does not have such property, (see Fig. [T]). 

Let us hrst briehy explain how the extended DEF model can be consistently embedded in the 
cosmology before describing quantitative analysis in the subsequent sections. As is the case with 
the original DEF model, in the extended DEF model the scalar held at the origin in the presence of 
matter becomes unstable and should in principle be pushed away from the origin. Thanks to the 
mass term, there exists a global minimum of the effective potential, which helps the asymmetron 
to settle down at this point. Assuming the universal conformal coupling, the scalarization should 
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happen during inflation. Due to the non-vanishing value of the asymmetron, the gravitational 
constant would be different from the one we measure in the laboratory and in this sense the 
law of gravity would be different from GR as we know. After inflation, the Universe is reheated 
and dominated by radiation. Since the trace of the radiation energy-momentum tensor is zero, 
asymmetron decouples from the matter and the global minimum shifts back to the origin of the 
effective potential. As the Universe further expands, the Hubble parameter gradually decreases 
and eventually becomes smaller than the mass of the asymmetron. By then, the asymmetron 
undergoes damped oscillation and the Universe gradually approaches GR. That is, GR is a 
cosmological attractor in this model. As a result, GR is recovered to a good approximation in 
the present Universe. We will further show that osillating component of the asymmetron, which 
interacts only gravitationally with standard model particles, is a good candidate for cold dark 
matter. Therefore not only is our extended DEF model cosmologically viable but also it provides 
a dynamical mechanism for dark matter generation via asymmetron production during inflation. 

Of course, cosmology is not the only arena where the extended DEF model becomes relevant to 
observations. When the matter density inside a compact astrophysical object exceeds the critical 
density, asymmetron would undergo the spontaneous scalarization and the laws of gravity may 
deviate from GR by a considerable amount. This phenomenon itself is similar to the original DEF 
model, but our asymmetron model provides additional new features as follows. Due to the mass 
term, asymmetron outside the compact object where no matter exists diminishes exponentially 
on the length scale of the inverse of the mass. This is in clear contrast to the original DEF model 
where the asymptotic value is arbitrary and fixed by the boundary condition. Furthermore, the 
critical density beyond which the spontaneous scalarization occurs is not necessarily around the 
matter density of the neutron star (it could be either higher or lower), whereas in the original 
DEF model the spontaneous scalarization occurs inside the neutron star where the gravitational 
energy becomes comparable to its rest mass energy. This opens up a new possibility that not 
only neutron stars but also less compact astrophysical objects are the best targets to search for 
the deviation from GR. 

In this paper, we analyze in detail the spontaneous scalarization in the asymmetron model 
and how the gravitational field changes outside the compact star before and after the sponta¬ 
neous scalarization. We also show that inflation, assuming the universal conformal coupling to 
all the matter fields, induces the spontaneous scalarization and the asymmetron undergoes co¬ 
herent oscillations in later time Universe. As mentioned above, such an oscillating field can be a 
candidate for cold dark matter. We show that there is a parameter space where the production 
of asymmetron can saturate the dark matter content. In the last section, we further comment 
on the possibility of the asymmetron as dark energy. 

2 Spontaneous scalarization in high density region 

2.1 Model 

We introduce a real massive scalar field (p which is universally coupled to all the matter fields 
including the standard model particles through the metric g^i, = A^{4))g^y (Thus, g^j_u is the 
Jordan metric). This ensures that weak equivalence principle is satisfied. We assume that g^^u 
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satisfies the Einstein equations. Therefore, the basic action is given by 


S — + S^[g, (j)] + Smidfj.i-'i 4^] 

= j d'^Xy/^ ( iGTrGjv ~ ~ + / d‘^Xy4^C„,{g,'llJm), (1) 

where Gn is the Newton’s constant and Cm is the matter Lagrangian of all the matter fields 
including the standard model fields. The corresponding equations of motion are given by 


□ (/) - + A^{(j))A^^f = 0, 


Gni/ — SttGn 


- ( ^9^''dg(t)dy(l) + ) g^u + dg(l)du4’ + A^{4>)Tgy 


( 2 ) 

(3) 


where 



2 6Sm 


(4) 


is the energy-momentum tensor with respect to g^u and T = g^'^T^y. Since all experiments are 
done with respect to g^^y, T^^y is the normal energy-momentum tensor we use in the standard 
general relativity. For this reason, we call T^y physical energy-momentum tensor 
independent of </>, the first EOM states that the effective potential for (p is 


#2 


Since T is 


V,si(P) = ^P^-'^A\p)f, (5) 

for which we have = 0. 

Now, let us choose the function A‘^{p) such that it is an even function and it monotonically 
decreases for p > 0 and asymptotically approaches a constant value. One simple form that 
satisfies all these conditions is given by 


A^P) 


I — e + ee 



( 6 ) 


with 0 < e < 1. Throughout this paper, we consider this form of A{p) and assume that e = 0(1) 
and is not very close to 0 nor 1. With this choice, the effective potential for p in the presence of 
non-relativistic matter, for which T = —p, becomes 

-e + £e"^^ p. (7) 

The shape of 14fj is shown in the right panel of Fig. [TJ When p is uniform and the system is 
static, p would take a constant value p which minimizes the effective potential. Taylor-expanding 
Ves{p) around p = 0, we have 

VesiP) = i ^2 ^ 0(^4)^ (8) 

it should be, it can be verified by explicit computation that the conservation law = 0 is an automatic 

consequence of the combination of Eqs. © and ®. 
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We find = 0 is stable for p < ppT = 2/r^M^/e, but becomes unstable when p exceeds ppT- 
When p > /OpT, ^ is given by 

^ = ■ (9) 

We see that (j) depends only logarithmically on p. Thus, unless p takes extremely huge values, 
(j) is 0{M). To conclude, the scalar field undergoes spontaneous scalarization when p > ppp is 
realized. 

Let us next consider how the gravity behaves in the symmetric phase where (/) = 0 is satisfied. 
In the symmetric phase, the interactions between (j) and the other matter fields are, in the leading 
order, written as ~ cjP'T^i,/M'^. As we will see later, for observationally interesting cases, M is 
typically much larger than TeV scale, he., far beyond the energy scale accessible by any terrestrial 
experiments. In this sense, (f) completely decouples from the other matter fields and behaves as 
a free massive scalar field. Since = 1, assuming there is no excitation of (j) field, Eqs. ([3j) 

reduce to the Einstein equations. If 4> is excited around </> = 0, excitation will be observed as 
dark component interacting only gravitationally with ordinary matter. It is then natural to 
suppose that such excitation constitutes (a part of) dark matter. More detailed analysis of this 
possibility including its production mechanism will be discussed later. Therefore, at low density 
region p < ppp, GR is recovered. 

2.2 Gravity in spontaneous scalarization phase 

Contrary to the symmetric phase, deviation from GR occurs in the scalarization phase, which we 
will investigate in the following. In the scalarization phase, due to a non-vanishing (j), matter fields 
interact with (j) with interaction strength proportional to (p. This acts as a fifth force between 
matter fields. Since p is massive, the interaction range of the fifth force is limited to lip. 
In addition to the emergence of the fifth force, field equations for gravity are also modified. 
Assuming no excitation of (p field around (p, Eqs. ([3]), rewritten in terms of the Jordan-frame 
metric become 

Gfxu T (fO) 

where 


f{^,PPT/p) 

o., = AH4')aN = (i - ppt/p) ) 

We find that satisfies the Einstein equations with the gravitational constant replaced by 
Geff and with the effective cosmological constant Agg. Contrary to the case of the standard 
Higgs mechanism, for which smaller cosmological constant is realized in the symmetry-breaking 
phase compared to that in the symmetric phase, opposite phenomenon happens in the current 
model. Namely, if there is no (or very tiny) cosmological constant in the symmetric phase, then 
a positive vacuum energy of O(ppt) emerges in the spontaneous scalarization phase. Therefore, 
if the matter density is larger than ppx but is still the same order of magnitude of ppT, the 



Agff = AttGnP^P'^A ^{(p) = AttGnsppt In /(e, ppt/p) {l-e + 
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effective cosmological constant will also play a non-negligible role of gravitational physics. In the 
very high density region in which p S> ppT; we have 


Aeff ~ AttGn- - -p-pT In [ (1 — e)- ^ 


1 — e 


Pvt 


Geff ~ (1 — s)Gn- 


(13) 


We find that Aeg is enhanced only logarithmically from ppT- Thus, in the very high density 
region, effect of the effective cosmological constant is much smaller than the right hand side 
of (jlOp and does not significantly affect the dynamics. The effective gravitational constant is 
reduced by (1 — s). Thus, gravity is weakened by this amount. 

In the above argument, we have ignored the contribution of the scalar force and focused only 
on the change in the pure gravity sector. In order to evaluate the scalar force, let us consider a 
test point source of its physical mass Ms immersed in the static and uniform matter distribution 
in which spontaneous scalarization occurs. Presence of the point source distorts the scalar held 
from 0 by the amount 6(j) as well as the Einstein-frame metric from by the amount h^i, 1^^ . 
We assume Ms is so small that both 5(f) and can be obtained by linear perturbation analysis. 
Then, equation for 5(f) is obtained by linearizing Eq. ([2]) on the background (f> = (p given by 
Eq. Q. On this background, we have 


A^AJ.7 = - 


A^^ 
M ’ 


where A = A{(p) and we have introduced a dimensionless parameter ^ dehned by 

. ^ £ Vln/(e,ppT/p) 

~ \/2 f{£,pPT/p) 


(14) 


(15) 


Notice that in the deep scalarization phase for which p 3> ppp, this parameter is suppressed by 
a small factor pptIp- Since ^ controls the conpling between the asymmetron and matter helds, 
the coupling is weak in the deep scalarization phase. 

Using this quantity, the equation for 5(f> becomes 


(A - p^) 5 (f) = - 


-jfPs, PS = 


(16) 


Solution of this equation is given by 


5(f){r) 


dvr M Ar 


(17) 


Metric perturbation can be obtained in the standard manner. Noting that the gravitational 
constant is A'^G^ in the scalarization phase and Ar is the physical distance, we have 


/iQo — 217, hij — 21/5, 




U = 


AGnMs 


in the isotropic coordinates (or the PPN coordinates) 


#4 


(18) 


For simplicity, we do not take into account the cosmological constant term given by Eq. El) which exists in 
the background. Inclusion of it is straightforward. 

’^^If we are living in the scalarization phase, we have to replace MGn by Gjv to satisfy the local gravity 
experiments. See the last paragraph of the last section for relevant discussion. 
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The Jordan-frame metric with first order deviation from the background is given by 

dtJ-u — {4’)9tJ,u = A + h^i, -\- (Inyl ) ^1^=^ ■ (19) 

Since the constant overall factor A'^ is irrelevant to the following discussion, we will omit it. 
Substituting the above results to we find 


900 — — 1 + 2 [/ -|- 


2'kM'^ r 


9ij 


(^1 + 2U 


2'kM‘^ r ) 


5ij. 


( 20 ) 

( 21 ) 


We find that the scalar force described by the Yukawa potential contributes to the metric per¬ 
turbation in the Jordan-frame which does not match the form predicted by the pure GR. We can 
translate this contribution to the PPN parameter 7 (see, for instance, [T]). This parameter is 
defined by gij = {l + 2'^JJ)5ij where IJ is metric perturbation of the 00 component, 500 = —1 + 2[/ 
(in GR, 7 = 1 ). In the present case, 7 becomes 


2A ^ 

1 + A’ ^ ' 


( 22 ) 


Since M appears in the denominator of in Eq. ([2]), naively one would expect that if M 
is comparable or smaller than the Planck scale ~ ' , then the scalar force would become 

stronger than the gravitational force within the range ~ The above result shows that this 
naive expectation is not correct since it is ^^/(M^Gat) that determines the magnitude of the 
deviation from GR. As we mentioned earlier, ^ becomes small in the deep scalarization phase 
and the system can become close to GR (I 7 — 1| <C 1) even when M < G^ ' . 


2.3 Spontaneous scalarization only inside a compact object 

Having explained the basic picture of the spontaneous scalarization, it is intriguing to analyze a 
situation where a dense object inside which spontaneous scalarization occurs resides in vacuum. 
To capture the essence of the phenomena, we make the following simplification that the object 
is static, uniform and spherically symmetric and is made of non-relativistic matter and its size is 
much larger than the Schwarzschild radius so that metric in the Einstein-frame can be taken to be 
the Minkowski one, but density is much larger than ppT- These assumptions will be inappropriate 
in quantitative sense for dealing with realistic astrophysical objects such as normal stars, white 
dwarfs, neutron stars and so on. But we believe that the following result remains qualitatively 
correct. 

With the above assumptions, the equation for (j) becomes 

^ r dr d4) ^ ^ 

As is done in [36], let us perform the change of variables as 

r ^ r, 4>^ X, Ves ^ -U. (24) 
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Figure 2: Potential ?7 for r < i? and t > R. Initially, x is at Xc and asymptotically approaches 
zero for r —>■ oo. 


Then, the above equations becomes 


(fx 2 dx 


dU 


dr'^ r dr dx' 


(25) 


which represents motion of a point mass under the potential U associated with time dependent 
friction. Denoting R by the radius of the object, U changes its shape at r = i? as shown in 
Fig. [2j What we want is a solution x{t) with boundary condition. 


a:(0) = Xci i;(0) = 0, x{t —>■ oo) = 0. 


(26) 


We follow [36] to construct the approximate analytic solution for this kind of problem. 

When R is large enough, x stays near Xc for a long time. This means that Xc is very close 
to (p at which U' = 0. At t = R, the friction had become negligible and the kinetic energy of 
X is just enough to be compensated by the difference of potential energy between x = x{R) and 
X = 0 so that X asymptotically approaches zero. Since the transition at t = R happens near Xc, 
let us replace U (x) before t = Rhy quadratic form around 4>, 


U{x) « -^m‘^{x - (pf, 

where m? is the second derivative of I4ff evaluated at cp, 

2/^2 


— Ves,(p<j){(p) = —^ ( 2ppT — (1 — s) 


PPT V 


/(Ej^pt/p) 


ln/(e,/9pT/p). 


Then, Eq. (1251) becomes 


2 dx n, 

TT + = m [x- (p), 

dr^ T dr 


(27) 


(28) 


(29) 


whose solution with the initial condition x(0) = x^-, ic(0) = 0 is given by 

/ \ 7 / sinh(mT) 

x{t) = (p+{xc- (p) -• 


mr 


(30) 










On the other hand, x{t) for t > R with boundary condition x —>■ 0 for r —>■ oo is given by 


x{t) = C - 


(31) 


where C is integration constant. Requiring x and x are continuous at t = R determines Xc and 
C as 


—m(l + R^) + m cosh(mi?) + ^ sinh(mR) - 

mcosh(mR) + /isinh(mR) 

mi?cosh(mi?) — sinh(mi?)- 

C = --r-- rCp. 

mcosh[mR) + fismh[mR) 

In the high density limit p ^ Ppt, becomes 


= 2p^ In 


(1 

PPT 


(32) 

(33) 


(34) 


which is enhanced by log factor compared to /i^. Then, neglecting p terms in Xc yields 

-(1 + Ru) + cosh(mR) - 
cosh(mrij 

Thus, if i? ;§> p~^) is satisfied, then cj) stays very close to ^ until the surface of the 

object and then decays exponentially over length scale p~^ outside the object. In other words, 
we can say that spontaneous scalarization occurs inside the object when the size of the object is 
much greater than the Compton wavelength of </> in symmetric phase (in addition to the trivial 
condition that density is greater than ppt)- 


(35) 


2.4 Gravity outside the scalarized compact object 


Let us consider metric perturbation outside a compact object inside which spontaneous scalar¬ 
ization occurs. As in the previous subsection, we assume that the compact object is made of 
non-relativistic matter. We assume that matter density is high enough so that spontaneous 
scalarization occurs inside the object but not compact enough so that gravity is weak every¬ 
where. From Eqs. jSl), we see that this amounts to perform perturbative expansion of the metric 
in the Einstein-frame around the Minkowski metric in terms of a dimensionless quantity given by 
(Schwarzshild radius)/(distance) 1^. In this subsection, we consider only first order correction 


and treat the linearized Einstein equations. 

One may wonder why we do not consider linear perturbation in the Jordan-frame. In order to see this is not 
feasible, let us express the field equations m in terms of the Jordan-frame metric. They are given by 


InTVd In A—2V“Vci lnA)+2\7^Vi, In A = 8ttGn 




We find that terms containing In A on the left hand side of the above equation are not associated with Gn- This 
makes sense since they come from for the Einstein-frame. It is now clear that Jordan-frame metric cannot be 
expanded in terms of (Schwarzshild radius)/(distance). Indeed, since A changes by 0{e) from inside to outside of 
the compact star, a large variation of the Jordan-frame metric (exceeding CJ(Gjv)) is induced near the surface of 
the compact object. 
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We decompose the metric in the Einstein-frame as 




(36) 


where hf^i, is proportional to Gat. As usual, we introduce by and impose 

a gauge condition d^h^y = 0. Then the linearized field equations become 


^ diag(A"'p, 0, 0, 0). 


= -167rGArA^(()))T),j,. 

For the non-relativistic matter, we have 

T, 

Thus, only the t — t component becomes non-trivial, 

Ahoo = — 

Then, the gravitational potential defined by = — ^/iqo becomes 


<h(x) = —Gat / ux' 


/ 


,3 , A^p{x') 


\x — x'\ 

In particular, when x is very far from the object, this becomes 


‘h(x) « —9IL j (fix' A^p{x'). 


(37) 

(38) 

(39) 

(40) 

(41) 


Since symmetry is restored {4> = 0) outside the star, Einstein-frame is equivalent to the Jordan- 
frame in such a region. Thus, the physical gravitational potential $ is also given by Eq. m- 
The distance r approaches the physical distance when r is much larger than the size of the star. 
Noting that A^d?x is the physical volume element, physical mass Ms is given by 


Ms = / d^x A^p{x). 


(42) 


If the size of the star is much bigger than spontaneous scalarization occurs inside the star 
and (j) takes the uniform value (p given by Eq. ([9]) everywhere inside the star except for the thin 
shell region near the surface. Thus, it is reasonable to approximate A to be uniform inside the 
star {A = A{(f>) = Am) and to have a step-function like transition at the surface of the star and 
to become unity outside the star. With this simplification, we have 


‘h(x) « — 


GnAi^Ms 


(43) 


Eor p ppT, we have Am ~ — s. Therefore, from the observer outside the star, M appears 

to be decreased by Ain, or equivalently, Gn appears to be decreased by Ajn. 

Taking a component parallel to of the conservation law = 0 for the non-relativistic 

matter and for the metric = A^r/^j^, we have 




(44) 


10 





where is defined by = v'^/A. Thus, the mass Ms defined by Eq. (|42l) . which is the sum 
of mass of each particle that constitutes the star, is conserved unless no matter escapes/enters 
the star. This implies that gravitational potential far from the star changes by after the 
star undergoes the spontaneous scalarization. At first glance, it appears that this conclusion 
is inconsistent with the Birkhoff’s theorem. In order to understand this in more detail, let us 
consider a spherically symmetric star whose density is initially smaller than p-pT- Let us assume 
that, at some time for some reason such as reduction of the radiation pressure due to depletion 
of fuel to produce thermal energy, the star starts to shrink and the density eventually exceeds 
Ppt before the star settles down to a new stable configuration. By the time when the star 
becomes static again, the spontaneous scalarization is realized inside the star. This final state is 
already described in previous subsections. Let us write the Einstein-frame metric describing the 
transition as 

ds^ = —(1 + 2<I>(t, r))dt^ + (1 + 2A(t, r))dr‘^ + (45) 

where both 4> and A are treated as linear perturbations just as the previous subsection. The 
scalar field also respects the spherical symmetry and hence (j) = (p{t,r). Outside the star, the 
t — r component of the Einstein equations ([3|) become 

A = 47rGNr^(j)'. (46) 


By integrating this equation along time with fixed r, we have 
r [A(f oo, r) — A{t ——oo, r)] = dvrGArr^ 


dt ^(j). 


(47) 


Erom the argument of the previous subsection, the left hand side of the above equation is equal 
to (Ain — l)Gj^Ms when r is much bigger than the radius of the star. Thus, we have 


(Ain - l)Ms = dtSr 


(48) 


where Sr = 47rr^ is the surface area of the sphere of radius r. This result shows that change 
of the gravitational potential before and after the spontaneous scalarization is compensated by 
the emission of the scalar wave whose flux is given by Neglecting the metric perturbation, 
equation of motion for (j) is given by 


— (j) + 4>' — 0 . 

r 


(49) 


Before the star starts to shrink, since the density of the star is less than ppT, </> = 0 everywhere. 
After the star starts to shrink and when the density exceeds ppT, at (^ = 0 becomes unstable 
inside the star and this acts as a force to push (p into the stable point. In this way, p inside 
the star changes its value 1^. This change also excites the change of p outside the star and it 


Since scalarization with positive cj> and negative one are equally allowed, scalarization occurs randomly on 
distance over the correlation length. As a result, the compact star just after the spontaneous scalarization may 
be mixture of positive and negative p and two regions are separated by domain wall. Though this may lead to 
interesting phenomena, process of spontaneous scalarization with this effect being taken into account is complicated 
and we do not consider it in this paper. 
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propagates as a wave which decays as ~ 1/r. Contribution of the scalar wave to the Jordan-frame 
metric far from the star is given by 




-2 


(50) 


Thus, this contribution is more suppressed for large r compared to gravitational potential and 
gravitational wave both of which decay as ~ 1/r although the latter is absent in the present case 
from the beginning due to the simplified assumption that the system is spherically symmetric. 
For a distant observer, the dominant deviation from GR caused by the spontaneous scalarization 
is the change of the gravitational constant. 


3 Asymmetron as dark matter 


Having introduced a new scalar field (j) which interacts with standard matter only gravitationally 
in the symmetric phase, it is natural to identify it with dark matter. As we will show, the 
spontaneous scalarization also provides a natural mechanism of fixing the abundance of dark 
matter within the framework of primordial inflation. 

Let us consider the effective potential during inflation. Making the phenomenological ap¬ 
proximation that inflation is caused by the fluid with its equation of state i/nf = —pini, we 
have ^ 

+ Pinf. (51) 

True effective potential differs from this potential by the amount of slow-roll parameters multi¬ 
plied to the second term, which is small enough for our present purpose and we ignore it. Due 
to the contribution of the pressure, the coefficient of the second term on the right hand side is 
enhanced by a factor four compared to the case of the non-relativistic matter. As a result, (p 
when the spontaneous scalarization occurs is given by 

= ln/(e,ppT/(4pinf)), (52) 


From this, we find that spontaneous scalarization occurs for > ppt/J, which we assume to 
be satisfied. The critical density is not equal to ppT because ppT is defined as the critical density 
for the case of the non-relativistic matter (see below Eq. ([8|l). 

During inflation, the (p held is hxed to the value given by Eq. (I52p . On top of this, since (p 


is almost massless during inhation Eeff,(/>(/> ^ perturbations originating from the quantum 


huctuations of (p are generated. We will come back to this issue later. After inhation, the Uni¬ 
verse is reheated and dominated by radiation. When this happens, T vanishes and the effective 
potential reduces to the bare potential 1^1 . As the Universe expands, the Hubble parameter 


’^^Strictly speaking, this is not correct since there exists non-relativistic baryon component even in radiation 
dominated era after the QCD phase transition which occurs around temperature Tqcd ~ 200 MeV. The baryon 
density at this temperature is estimated as pb(TQCD) ~ 6x 10“^^ GeV"^ for fit, = 0.04, g^* = 20, Tqcd = 200 MeV. 
If ppT is smaller than pbiTQco), which is the case for p < 10“^^ eV (see Fig. [3]) when we require asymmetron to 
be dark matter, the baryon forces the asymmetron to undergo the spontaneous scalarization at the time of the 
QCD phase transition. As a result, the result (I56II cannot be applied straightforwardly and we need to modify it 
in an appropriate way. Since pbfTqcD) is much smaller than the nuclear density which is an interesting target for 
ppT, we do not consider this case in this paper and set p > 10~^^ eV. 
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gradually decreases and at some point becomes equal to /x. Before this time, cf) keeps its initial 
value fixed during inflation. After this time, (j) oscillates around the origin like ~ sm{pt)/t 
and behaves as non-relativistic matter. Thus p^ decreases as 1/a? in the Einstein-frame (a is 
the scale factor in the Einstein-frame). Then, the energy density of (p at present time is given by 


Ppo 


"/^0,OSC 


^eq 


(1 + ^^eq) 


1 /y^ — 

_ _ _ ^osc M i 2 

(l + Zeqf a3q 2*^ 


(53) 


where the subscript osc in any quantity means that it is evaluated when p starts oscillations, 
i.e., p = Hose and Ueq = 1/(1 -|- Zeq) is the scale factor at the time of matter radiation equality. 
We assume that there is no additional entropy production after inflation. Therefore, the entropy 
density of radiation decays as 1/a? = A^a? (a is the scale factor in the Jordan-frame). With this 
assumption, we have 




5*s,eq / 9*,osc 


3/4 


( £l 

(1 J- ^^eq)^ 9*s,osc V 9*,eq J \Pr, 


eq 


3/4 


12 


2 


(54) 


where g*s, 9* represents the effective degrees of freedom entering in the entropy density, energy 
density of radiation, respectively. By the time of matter-radiation equality, amplitude of p has 
decreased enough so that there is little difference between the Einstein-frame and the Jordan- 
frame, i.e. ,Aeq ~ 1 to a very good approximation. Now, by using the Friedmann equation in 
the Einstein-frame 


..2 _ tj2 _ SttGtv ^4 ~ 

P ~ ^ osc ~ Q ^infPr,osc, 


(55) 


to eliminate /5r,osc and Eq. (j52p to eliminate p, we end up with 




P(l>,0 _ 9*s,eq I 9*,osc 


Pc,0 9*s,osc \ 9*,eq 


3/4 


eppT 


In / (e, ppT/(4pinf)) i ( — 


^inf 


TT \ 3/2 
Ho \ 3/4 

P J 


(56) 


If this quantity is equal to the observed then p constitutes the whole dark matter. This 

requirement yields a relation between p and ppT, which is shown as ^dm = Mdm(ppt) in Fig- El 
In this figure, the parameters are hxed as 5'*^eq = 9*s,eq = 100, and e = 1/2. 

Since there is strong upper limit on the deviation from general relativity by the solar system 
experiments as well as terrestrial ones, we require that the spontaneous scalarization occurs at 
density larger than the Earth density, i.e., ppx PEarth ~ 5 x 10“^^ GeV^Combining this 
with the footnote #7 our primary interest for ppp is ppT ^ 10~^^GeV^. Then, from Fig. [3l we 
find that the corresponding restriction on p is given hy p > 10“^^ eV which we regard as the 
possible minimum value of our interest. 


3.1 Isocurvature constraint 

We saw in the previous subsection that spontaneous scalarization occurs during the primordial 
inflation and this provides a mechanism of preparing non-zero value of the asymmetron to realize 

#®It is possible that ppx < psarth and we are living in the spontaneous scalarization phase. One possibility is 
that e is very tiny e ^ 1. Since large deviation from GR never happens in any situation for such a case, we do 
not consider this possibility in this paper. The second possibility is that ppx is the order of the critical density 
of the Universe. In this case, asymmetron behaves not as dark matter but as dark energy. We will briefly discuss 
this scenario in the last section. 
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its coherent oscillations which behave as non-relativistic matter interacting only gravitationally 
with other matter fields. There is indeed a parameter range of /U and M where the energy density 
of asymmetron is equal to that of dark matter. But before asymmetron can be considered as 
a candidate of dark matter, it must satisfy other observational constraints. There are two non¬ 
trivial observational constraints, which we will consider below. 

The hrst constraint is the non-detection of the dark matter isocurvature perturbation. Since 
the (p field is almost massless during inflation, this field acquires almost scale-invariant classical 
fluctuations during inflation when each wavelength mode crosses the Hubble horizon. In addition 
to this, the standard adiabatic perturbations are also generated from classical fluctuations of 
either inflaton or other light helds, which are equally shared by all the existing particle species 
such as photons, baryons and dark matter. On top of this, dark matter has its own fluctuations 
coming from the fluctuations of the p field itself and these fluctuations act as isocurvature 
perturbations having no correlation with the adiabatic ones. Since there is a strong upper limit 
on the amplitude of the isocurvature perturbations imposed by CMB observations, this limit 
can be converted to the constraint on the domain of the (/U, /Opt) plane. To see this in more 
quantitative manner, let us first introduce the dark matter isocurvature perturbation 5 dm by 137! 


5dm 


i^Pdm _ 3 6 p ^ 
Pdm 4 


(57) 


where 6p^ is the density perturbation of photons. This quantity is conserved as long as the scale 
considered is super-Hubble scale. In the present case, /Odm = Pcf)- 

As is done in [38] (but for a different model), we adopt the approximation of the sudden 
transition that p is completely frozen before H = p (hence P(j, is constant in time) and starts 
to oscillate exactly when H = p and behaves as non-relativistic matter p^ oc [39]. Then, 
the hypersurface on which p starts to oscillate coincides with the one with constant total energy 
density. Since 5dm is independent of the choice of time slicing, we can compute 6r and 5^ in any 
time slicing and we take the H = m hypersurface for this purpose. On this hypersurface, we 
have 

3^2 

Pr{x) + pcpix) = ptot = ( 58 ) 

Decomposing this relation into the background part and perturbation part and extracting the 
perturbation part, we have 

<5r(^) = (59) 

where = Pcf,/ptot, evaluated at time when H = p. Plugging this relation into the definition of 

5dm ) we have 

Sdm{x) = ~ (®®) 

where we have used ^ 1 since the time when asymmetron starts to oscillate for the range of 
p of our interest is much earlier than the time of matter-radiation equality. Since p^ = pP(p?‘ j2 
in the radiation dominated era, we finally have 


5dm (^) 


2bp 

P 


(61) 
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Here, 8(j) is the perturbation quantum mechanically generated during inflation. This is un¬ 
correlated with the (adiabatic) curvature perturbation which is sourced by other fields. The 
corresponding power spectrum of 5 dm is given by 

— — ( ^ ^ Pinf 


where the modified Friedmann equation 


1 


Svr 


^inf 


e/3PT In /(e, ppx/(4/iinf)) ’ 


= 87rGArA'^((())Anf, 


(62) 


(63) 


is used to obtain the final expression. The upper limit on the uncorrelated dark matter isocur¬ 
vature perturbation by WMAP 9yr is given by m 

a < 0.047 (95% CL), (64) 

Vn I - a 

where V-r, is the power spectrum of the adiabatic perturbation. For fixed ppp and pinf, this 
bound can be converted into the upper bound on // which is shown as a line of p = piso in Fig. [3l 
We find that a line p = pdm lies above p = piso for the case of high energy inflation scale such 
as pinf = (10^^ GeV) . Thus, the scenario of asymmetron being responsible for the total content 
of dark matter is incompatible with those inflation models where the energy scale is as high as 
~ 10^^ GeV. The line p = pdm comes below the isocurvature constraint line if the inflation 
energy scale is lowered, which can be understood from the expression of Pdm given by Eq. (1621) . 
The equation shows that Pdm is basically proportional to pinf (the denominator depends only 
logarithmically on Pinf). Indeed, if pinf is as low as 10^^ GeV, the isocurvature constraint is 
evaded for all the range of ppp we are interested in {i.e. nuclear energy density). 


3.2 Constraint from the fifth force experiments 

As is already mentioned, we are interested in the case where ppp is between the stellar density 
and the nuclear density realized in the neutron stars so that spontaneous scalarization takes 
place in compact astrophysical objects and 0{e) deviation from general relativity occurs only in 
such regions. In the asymmetron model, the solar system is in the symmetric phase (c^ = 0). As 
we saw in the previous section, GR is exactly recovered in this phase and this model passes the 
solar system and terrestrial experiments that have placed very strong limit on deviation from 
GR. However, this conclusion must be reconsidered more carefully if we take the scenario of the 
asymmetron being dark matter. In this case, cj) is coherently oscillating in time with angular 
frequency p, which describes the cold dark matter. The value of (j) averaged over time longer 
than the oscillation period is zero, but the value at each time is different from zero. Therefore, 
the assumption of no excitation of cj) in the symmetric phase is violated if we require cj) to be 
dark matter. 

When (j) is oscillating in the symmetric phase, the interaction between matter fields and (j) 
field given by A^A^T (see Eq. ([2])) also oscillates as —£asm{^t)/M‘^ (A ss 1 is assumed), where 
a is the amplitude of the oscillations of 4> and is given by ((•••) represents the 

time average over oscillation period). The amplitude a is determined from the requirement that 
= p^(</>^) coincides with the dark matter density. This condition yields 

2 2pDM 
a = — TT- 
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p,„fV4=ioi3 GeV 



(a) p\it = 10“ GeV 


p,„fV4=ioi= GeV 



(b) p\it = 10“ GeV 


Figure 3: Three curves in each panel are /tdM) l^iso and /xsth- Right to each curve of //iso and 
/isth is the allowed region satisfying the observational constraints. 
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The effect of oscillating cj) can be observed as the periodically time varying gravitational force 
(period is vr///) acting on two massive bodies with maximnm given by ~ ea/M^ on top of the 
standard gravitational force. In order to see this, let us determine the gravitational potential 
in the Jordan-frame. To simplify the analysis, based on the fact that the time scale of interest 
(e. g., time duration of experiments) is much larger than the oscillation period 2Tr/fi; 


27r 


4 X 10"^ s 




10-11 eV 


( 66 ) 


for a range of interest /r > 10“^^ eV, we make the approximation that only the averaged value of cj) 
enters the measurable gravitational potential and asymmetron both of which are simultaneously 
generated by the source object such as the Earth. Denoting 64> by the small deviation from the 
background (j) caused by the presence of the point mass with its mass Mg {p = Ms6{x)), equation 
for 54> is obtained by linearizing Eq. ([2]); 


Adcj) — p?5(f) + 


'_\/W 

2 M2 


-p = 0. 


Solution of this equation is given by 


(67) 


^ e^W) Ms 
^ 87rM2 


As a result, time-time component of the Jordan-frame metric is given by 


( 68 ) 


500 = -1 + 


2GjvMg 

r 


Fir), 


F{r) = I + 


^-ur 

WttM^Gn 


(69) 


The function F{r) represents the modification from the standard gravitational potential^. Elim¬ 
inating (())2) by Eq. ([65]), we have 


F{r) 


P^ PPM 

AtiGnS /9p.p 


(70) 


As expected, the (p field contributes to the force between two bodies as Yukawa type force. 
Various experiments have been performed to test the inverse square law of gravity. One of the 
typical modifications of the inverse square law which is actively tested by experiments is exactly 
the form given by E(r). In |40] . deviation from the inverse square law of a form F{r) = l + ae~'^/^ 
is assumed and upper limit on a is summarized for a wide range of A from 10“® m to 10^^ m. 
We converted this constraint in A — a plane to the constraint in ppx — p plane. The result is 
shown as a green(dotted) curve in Fig. |3l The region satisfying the fifth force experiments is 
right to the green(dotted) curve. We find that the constraint from the fifth force experiment is 
much weaker than the isocurvature one and is always satisfied for any interesting range of ppT- 

In addition to the modification by F{r), change of the Newton’s constant caused by the change of A {{A^) yf 1) 
also modifies the gravitational. However, this correction is negligibly small and we do not take this effect into 
account in Eq. (PI) 
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4 Discussion and conclusion 


We have proposed the asymmetron model, a class of scalar-tensor theories, in which the significant 
deviation from GR occurs only in high matter density region. This is an extended version of 
the Damour-Esposito model proposed in [18] by adding mass term and allowing the energy scale 
appearing in the conformal factor to differ from the Planck scale. We have shown that the 
asymmetron model can be consistently embedded in the cosmological framework. In particular, 
spontaneous scalarization caused by inflaton in a dynamical way provides the initial condition 
for the subsequent coherent oscillations of the asymmetron. The damped oscillation has nice 
properties in that it not only makes the asymmetron behave as cold dark matter but also makes 
GR a cosmological attractor. Oscillating asymmetron yields periodically varying fifth force but 
its magnitude is far below the current experimental sensitivities and the model we studied in 
this paper is in practice indistinguishable from GR in the present Universe except inside dense 
compact objects and easily passes the solar-system and terrestrial experiments. There is a range 
of parameter space where the asymmetron can saturate the whole dark matter component and 
at the same time significant deviation from GR in the present Universe occurs only inside the 
dense compact objects such as neutron stars. 

In the spontaneous scalarization phase, the gravitational constant becomes smaller than that 
in the symmetric phase, namely the value determined in laboratories. Thus, the gravity is 
weakened only inside dense compact objects, which is a dominant modification from GR. The 
scalar force also operates among matter with strength given by ~ (see Sec. 12.2p 

compared with gravitational force. In the deep scalarization phase where <C 1, the scalar 
force can become tiny even when M is smaller than the Planck mass Gj^ ' . Furthermore, the 
interaction range is limited by ~ and for the range of our interest this scale is rather short. 
As a result, weakening of gravity is the dominant feature representing deviation from GR when 
the density of compact object is much bigger than the critical density above which spontaneous 

scalarization occurs. This suggests that the size of the compact star in the asymmetron model 

— 3/2 

becomes larger than that in GR. Since the Ghandrasekhar mass is proportional to ' , we 
expect that the Chandrasekhar mass in our model should be larger than that in GR for compact 
stars undergoing spontaneous scalarization. 

There are many issues that we did not consider in this paper and deserve further investiga¬ 
tions. In this paper, we mainly focused on the mechanism of the spontaneous scalarization in 
the asymmetron model, its basic properties and embedding it in the cosmological framework. 
Obviously, the next thing to do is to investigate how to test this model in astrophysics, in partic¬ 
ular in connection to gravitational wave observations. In this context, it is first of all interesting 
to clarify how the stellar structure (such as mass-radius relation and the Chandrasekhar mass) 
in the asymmetron model is modified from GR. Gravitational waves from compact binaries are 
the main target of the laser interferometers. Dynamics of binaries, waveform of the gravitational 
waves and detectability by using the interferometers for the original DEF model have been stud¬ 
ied |241t33] . Performing the similar analysis for the asymmetron model will help to elucidate 
what observation is the best probe for exploring the asymmetron model. 

Another intriguing thought is the possibility of asymmetron being responsible for dark energy. 
In this paper, we have focused on the case where ppT S> PEarth and the spontaneous scalarization 
occurs only in extremely high density region. On the other hand, if ppx is the order of the 
current critical density of the Universe, we expect that the scalarization persists until present 
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epoch and the mass term approximately plays the role of the cosmological constant. 

Taking M to be Planck mass, this is achieved if ^ is chosen to be around the Hubble constant 
Hq. This means that locally, such as in the solar system, asymmetron mediates a long range 
force in addition to the gravitational force. However, since the solar system is in the deep 
scalarization phase, the scalar force is suppressed by the factor ^ given by Eq. m- Indeed, 
if we take ppT to be the present critical density of the Universe and p the density of solar 
wind (we assume one proton per cubic centimeter) and e = 1/2, we have ^ ~ 5 x 10“^ and 
1 — 7 ~ 2 X This value is much below the current constraint I7 — 1| < 10“® obtained 

from the time delay measurement [T]. Thus, the asymmetron as dark energy can safely pass 
the solar system constraints. The most characteristic feature would be time-dependence of the 
gravitational constant (see Eq. (fT^ i through the time-dependence of the matter density p due to 
the cosmic expansion. On cosmological scales, the gravitational constant gradually increases as 
the Universe expands and it is interesting to investigate how the large scale structure is affected 
by such a time varying gravitational constant. 
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